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We demonstrate entanglement distribution between two remote quantum nodes located 3 meters 
apart [J]. This distribution involves the asynchronous preparation of two pairs of atomic memories 
and the coherent mapping of stored atomic states into light fields in an effective state of near 
maximum polarization entanglement. Entanglement is verified by way of the measured violation 
of a Bell inequality, and can be used for communication protocols such as quantum cryptography. 
The demonstrated quantum nodes and channels can be used as segments of a quantum repeater, 
providing an essential tool for robust long-distance quantum communication. 



In quantum information scienee , distribution of en- 
tanglement over quantum networks is a critical require- 
ment , including for metrology [3|, quantum computa- 
tion 0, Q and communication 0,1^ ■ Quantum networks 
are composed of quantum nodes for processing and stor- 
ing quantum states, and quantum channels that link the 
nodes. Significant advances have been made with diverse 
systems towards the realization of such networks, includ- 
ing ions 0] , single trapped atoms in free space [1, Q and 
in cavities [Toj, and atomic ensembles in the regime of 
continuous variables 

An approach of particular importance has been the 
seminal work of Duan, Lukin, Cirac, and ZoUer (DLCZ) 
for the realization of quantum networks based upon en- 
tanglement between single photons and collective exci- 
tations in atomic ensembles [T2j . Critical experimental 
capabilities have been achieved, beginning with the gen- 
eration of nonclassical fields with controlled wave- 
forms [isj and extending to the creation and retrieval of 
single collective excitations [l^, [3 with high effi- 
ciency [19I [20| . Heralded entanglement with quantum 
memory, which is the cornerstone of networks with effi- 
cient scaling, was achieved between two ensembles [2l| . 
More recently, conditional control of the quantum states 
of a single ensemble |23,, 24] and of two distant en- 
sembles [2^ has also been implemented, as are likewise 
required for the scalability of quantum networks based 
upon probabilistic protocols. 

Our interest is to develop the physical resources 
that enable quantum repeaters [6], thereby allowing 
entanglement-based quantum communication tasks over 
quantum networks on distance scales much larger than 
set by the attenuation length of optical fibers, includ- 
ing quantum cryptography [2^ . For this purpose, her- 
alded number state entanglement [2l[ between two re- 
mote atomic ensembles is not directly applicable. In- 



stead, DLCZ proposed to use pairs of ensembles (C/^, D^) 
at each quantum node i, with the sets of ensembles 
{Ui}, {Z?i| separately linked in parallel chains across the 
network T?| . Relative to the state of the art in Ref. [2l[ , 
the DLCZ protocol requires the capability for the inde- 
pendent control of pairs of entangled ensembles between 
two nodes. 

In our experiment, we have created, addressed, and 
controlled pairs of atomic ensembles at each of two quan- 
tum nodes, thereby demonstrating entanglement distri- 
bution in a form suitable both for quantum network ar- 
chitectures and for entanglement-based quantum com- 
munication schemes. Specifically, two pairs of remote 
ensembles at two nodes are each prepared in an entangled 
state 21 1, in a heralded and asynchronous fashion [25| . 
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thanks to the conditional control of the quantum mem- 
ories. After a signal heralding that the two chains are 
prepared in the desired state, the states of the ensembles 
are coherently transferred to propagating fields locally at 
the two nodes. The fields are arranged such that they ef- 
fectively contain two photons, one at each node, whose 
polarizations are entangled. The entanglement between 
the two nodes is verified by the violation of a Bell in- 
equality. The effective polarization entangled state, cre- 
ated with favorable scaling behavior, is thereby compat- 
ible with entanglement-based quantum communication 
protocols p^ . 

The architecture for our experiment is shown in Fig- 
ure 1. Each quantum node, L (left) and R (right), con- 
sists of two atomic ensembles, U (up) and D (down), or 
four ensembles altogether, namely {LU, LD) and {RU, 
RD), respectively. Each pair is first prepared in an en- 
tangled state, where one excitation is shared coherently, 
by using a pair of coherent weak write pulses to induce 
spontaneous Raman transitions \g) |e) \s), as illus- 
trated in the inset in Fig. 1. The Raman fields {Ilu,^ru) 
from {LU, RU) are combined at the 50-50 beamsplitter 
BSjj, with the resulting fields directed to single-photon 
detectors. A photoelectric detection event in either de- 
tector heralds that the two ensembles are prepared. The 
remote pair of D ensembles, {LD, RD), is prepared in 
an analogous fashion. 
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FIG. 1: Setup for distributing entanglement between two quantum nodes (L, K) separated by 3 meters. The inset shows the 
relevant atomic levels for the 6S1/2 — > 6P3/2 transition in atomic Cesium, as well as the associated light fields. The ensembles 
are initially prepared in \g). Weak write pulses then induce spontaneous Raman transitions \g) |e) |s), resulting in the 
emission of anti-Stokes fields (fields 1) near the |e) — > |s) transition along with the storage of collective excitations in the form 
of spin-fiips shared among the atoms [l^ . With this setup, a photo-detection event at either detector D\a or D\h (Die or 
-Did) heralds entanglement between the collective excitation in LU and RU (LD and _RD) [Ul (see text). BSu and BSd 
are implemented using two orthogonal polarizations in one fiber beam splitter, yielding excellent relative path stability. A 
heralding detection event triggers the control logic to gate off the light pulses going to the corresponding ensemble pair (f/ or 
D) by controlling the intensity modulators (/M). The atomic state is thus stored while waiting for the second ensemble pair 
to be prepared. After both pairs of ensembles f7, D are entangled, the control logic releases strong read pulses to map the 
states of the atoms to Stokes fields 2 via |s) — > |e) \g). Fields 2lu and 2ld {2ru and 2rd) are combined with orthogonal 
polarizations on the polarizing beam splitter PBSl (PBSr) to yield field 2^ {2r). If only coincidences between fields 2l and 
2r are registered, the state is effectively equivalent to a polarization maximally entangled state. {\/2)l,r'- rotatable half-wave 
plates. 



Conditioned upon the preparation of both ensemble 
pairs {LU, LD) and {RU, RD), a set of read pulses 
are triggered to map the stored atomic excitations into 
propagating Stokes fields in well-defined spatial modes 
via \s) — > |e) — > \g) by way of a collective enhancement 



[T2| (see inset in Fig. 1). This generates a set of four 
fields denoted by {2lui'^ru) for ensembles {LU, RU) and 
{2ljj,2iijj) for ensembles {LD, RD). In the ideal case 
and neglecting higher-order terms, this mapping results 
in a quantum state for the fields 2 given by 



H2,u,2nu,2.o,2nn) = 3 {\^h^u\^h^u ± e"*- |1)2., |0)2h.) ^ ® (|0) 2.. 1 1>2.. ± e"'" |1>2.. |0>2«.) ^ ■ (1) 



Here, In)^; is the n-photon state for mode x, where 
X e {2lu,2ru,2ld,2rd}-, and rju {rjo) is the relative 
phase resulting from the writing and reading processes 



for the U {D) pair of ensembles [2l| . The ± signs for the 
conditional states U, D result from the unitarity of the 
transformation by the beamsplitters {BSu, BSd)- The 
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extension of Eq. ([T]) to incorporate various nonidealities 
is given in the Appendix. 



Apart from an overall phase, the state 
|V'2i[/,2H!7.2i,D,2nr>) Can be rewritten as follows: 



2lo \vac)2R 



±(|0)2.,|l)2.„|0)2„,|l>2«„±e'(''--"-)|l>2,,|0)2,„|l>2„.|0>2,. 



(2) 



r 



where \vac)2i denotes |0)2i„ |0)2iD ■ If only coincidences 
between both nodes L, R are registered, the first two 



terms (i.e., with e 



') do not contribute. Hence, 



as noted by DLCZ, excluding such cases leads to an effec- 
tive density matrix equivalent to the one for a maximally 
entangled state of the form of the last term in Eq. ([2]). 
Significantly, the absolute phases riu and rm do not need 
to be independently stabilized. Only the relative phase 
rj = rju — r]D must be kept constant, leading to 1/2 unit 
of entanglement for two quantum bits (i.e., 1/2 ebit). 

The experimental demonstration of this architecture 
for implementing the DLCZ protocol relies critically on 
the ability to carry out efficient parallel preparation of 
the {LU, RU) and {LD, RD) ensemble pairs, as well as 
to stabilize the relative phase ry. The first requirement is 
achieved by the use of real-time control, as described in 
Ref. [2^ in a simpler case. As illustrated in Fig. 1, here 
we implement control logic that monitors the outputs of 
field 1 detectors. A detection event at either pair trig- 
gers electro-optic intensity modulators (/M) that gate off 
all laser pulses going to the corresponding pair of ensem- 
bles, thereby storing the associated state. Upon receipt 
of signals heralding that the two pairs of ensembles [LU , 
RU) and {LD, RD) have both been independently pre- 
pared, the control logic triggers the retrieval of the stored 
states by simultaneously sending a strong read pulse into 
each of the four ensembles. Relative to the case where no 
logic is implemented, a 19-fold enhancement is obtained 
in the probability to generate this overall state from the 
four ensembles (see Appendix). 

The second requirement, for stability of the relative 
phase 77, could be accomplished by active stabilization of 
each individual phase ri(j,ri]j, as in Ref. [2l|. Instead of 
implementing this challenging technical task (which ul- 
timately would have to be extended across longer chains 
of ensembles), our setup exploits the passive stability 
between two independent polarizations propagating in a 
sing le interferometer to prepare the two ensemble pairs 
[27| . No active phase stabilization is thus required. In 
practice, we find that the passive stability of our system 
is sufficient for operation overnight without adjustment. 
Additionally, we implement a procedure that determin- 
istically sets the relative phase 77 to zero. 

We also extend the original DLCZ protocol by, as illus- 
trated in Fig. 1, combining fields {2lu, 2ld) and {2fiu, 



2rd) with orthogonal polarizations on polarizing beam 
splitters PBSl and PBSr to yield fields 2^ and 2^, re- 
spectively. The polarization encoding opens the possi- 
bility of performing additional en tang lement purification 
and thus superior scalability [H, [29| . In the ideal case, 
the resulting state is now effectively equivalent to a maxi- 
mally entangled state for the polarization of two photons 



l<.2je// (X \H2,)\V2R)±e^^\V2,)\H2R) 



(3) 



where \H) {\V)) stands for the state of a single photon 
with horizontal (vertical) polarization. The sign of the 
superposition in Eq. ^ is inherited from Eq. ([T]) and 
is determined by the particular pair of heralding signals 
recorded by(I?ia,-Dib) and {Die, Did). The entanglement 
in the polarization basis is well-suited for entanglement- 
based quantum cryptography [l^, , including security 
verification by way of the violation of a Bell inequality, 
as well as for quantum teleportation p^ . 

As a first step to investigate the joint states of 
the atomic ensembles, we record photoelectric counting 
events for the ensemble pairs {LU, RU) and {LD, RD) by 
setting the angles for the half-wave plates {^)l,r shown 
in Fig. 1 to 0° such that photons reaching detectors D26 
and D2d {D2a and I?2c) come only from the ensemble 
pair U {D). Conditioned upon detection events at Dia 
or Dih {Die or Did), we estimate the probability that 
each ensemble pair U, D contains only a single, shared 
excitation as compared to the probability for two excita- 
tions by way of the associated photoelectric statistics. In 
quantitative terms, we determine the ratio [2ll | 



PX.ll 



PX,10PX,01 



(4) 



where px.mn are the probabilities to register m photo- 
detection events in mode 2lx and n events in mode 2iix 
{X = {U, D}) conditioned on a detection event at Di. A 

necessary condition for the two ensembles {LX, RX) to 

(2) (2) 
be entangled is that h -^ < 1, where /i)^ — 1 corresponds 

to the case of independent (unentangled) coherent states 

for the two fields [2l|. Fig. 2 shows the measured /i^-* 

versus the duration tm that the state is stored before 

retrieval. For both U and D pairs, /i^^^ remains well 
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FIG. 2: Suppression ft'^' of the probabilities for each ensemble 
to emit two photons compared to the product of the proba- 
bilities that only one photon is emitted, as a function of the 
duration tm that the state is stored before retrieval. The solid 
line gives a fit for the U pair. Error bars indicate statistical 



visibility of V = 78 ± 3% in good agreement with the 
experimentally determined V = 75%. Finally, the fact 
that one of the fringes is inverted with respect to the 
other in Fig. 3B corresponds to the two possible signs in 
Eq. ([3]). As for 9l — 45° the measurement is sensitive 
to the square of the overlap ^ of photon wavepackets for 
fields 2u^D, we may infer S,u,d — 0.85 from the reduced 
fringe visibility {V = 55%) in Fig. 3B relative to A , if 
all the reduction is attributed to a nonideal overlap. An 
independent experiment for two-photon interference in 
this setup has shown an overlap ^ ~ 0.90 (see Appendix), 
which confirms that the reduction can be principally at- 
tributed to the non-ideal overlap. Other possible causes 
include imperfect phase alignment ^ and imbalance 
of the effective state coefficients (see Appendix). 

With the measurements from Figs. 2, 3 in hand, we 
verify entanglement unambiguously by way of the vio- 
lation of a Bell inequality ^3V\. For this purpose, we 
choose the canonical values, Ol — {0°,45°} and 9ji = 
{22.5°, -22.5°}, and construct the CHSH parameters 



below unity for storage times tm ^ 10 fis. For the U 
pair, the solid line in Fig. 2 provides a fit by the simple 
expression /i*^^) = 1 — Aexp (— (tm/t)^). The fit gives 
A = 0.94 ±0.01 and r = 22 ±2 fj,s, providing an estimate 
of a coherence time for our system. A principal cause 
for decoherence is an inhomogeneous broadening of the 
ground state levels by residual magnetic fields [30]. The 
characterization of the time dependence of /i^^^ consti- 
tutes an important benchmark of our system, as will be 
clarified shortly (also see Appendix). 

We next measure the correlation function E{9l,0r) 
defined by 

j^, r. n \ dac + C'bd ~ Cad ~ C'bc /_x 

E{0l,Or)^ . (5) 

Here, Cjk gives the rates of coincidences between detec- 
tors D2j and D2k for fields 2, where j, k € {a, b, c, d}, con- 
ditioned upon heralding events at detectors Dia, Du and 
Die, Did from fields 1. The angles of the two half- wave 
plates {^)l,r are set at 9^/2 and 9r/2, respectively. As 
stated before, the capability to store the state heralded 
in one pair of ensembles and then to wait for the other 
pair to be prepared significantly improves the various co- 
incidence rates Cjk by a factor that increases with the 
duration tm that a state can be preserved [l^ (also see 
Appendix) . 

Fig. 3 displays the correlation function i? as a function 
of 6r, with Si = 0° in A and 45° in B . Relative to Fig. 2, 
these data are taken with increased excitation probabil- 
ity (higher write power) to validate the phase stability 
of the system, which is evidently good. Moreover, these 
four-fold coincidence fringes in Fig. 3A provide a further 
verification that predominantly one excitation is shared 
between a pair of ensembles. The analysis provided in 
the Appendix with the measured cumulative ft.^^^ param- 
eter for this set of data, /i^^^ = 0.12± 0.02 , predicts a 



S+ = |£;(0°, 22.5°) + £:(0°, -22.5°) 

+£;(45°,-22.5°)-£;(45°,22.5°)| (6) 
S- = |£;(0°, 22.5°) + £:(0°, -22.5°) 

+£;(45°,22.5°) -£;(45°,-22.5°)| (7) 
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FIG. 3: Measured correlation function E{6l,0r) as a function 
of On with 6l fixed at A: 0° and B: 45°. The excitation 
probabilities for the ensembles are increased by ~ 1.5 times 
relative to Fig. 2, with each point taken for 30 minutes at 
typically 400/hour coincidence rate for each fringe. Error bars 
indicate statistical errors. 
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FIG. 4: Measured CHSH parameters S±, for the two possible effective states in Eq. (|3]), as functions of duration tm for which 
the first ensemble pair holds the prepared state. The excitation probabilities are kept low for high correlation (as in Fig. 2). A 
and B: binned data. The horizontal thick lines indicate the size of the bins used. C and D: cumulative data. The coincidence 
rate for these measurements is about 150/hour for each effective state. Error bars indicate statistical errors. 



for the two effective states 1-02 2 )e// ™ Eq. For lo- 
cal, realistic hidden variable theories, S± < 2 [3l|. Fig. 4 
shows the CHSH parameters S± as functions of the du- 
ration Tm up to which one pair of ensembles holds the 
prepared state, in the excitation regime of Fig. 2. As 
shown in the Appendix, the requirements for minimiza- 
tion of higher order terms are much more stringent in 
this experiment with four ensembles than with simpler 
configurations [25|. 

Panels A and B in Fig. 4 give the results for our mea- 
surements of S± with binned data. Each point corre- 
sponds to the violation obtained for states generated at 
Tm ± ^tm/2 {Atm marked by the thick horizontal lines 
in Fig. 4) . Strong violations are obtained for short mem- 
ory times, with for instance S+ = 2.55 ± 0.14 > 2 and 
S- = 2.61 ± 0.13 > 2 for the second bin, demonstrat- 



ing the presence of entanglement between fields 2^ and 
2fi. Therefore, these fields can be exploited to perform 
entanglement-based quantum communication protocols, 
such as quantum key distribution with, at minimum, se- 
curity against individual attacks [l3, [12] ■ 

As it can be seen in Fig. 4, the violation decreases 
with increasing tm- The decay is largely due to the time- 
varying behavior of /i^^^ (Fig- 2 and Appendix). In addi- 
tion to this decay, the 5*+ parameter exhibits modulation 
with Tm- We have explored different models for the time 
dependence of the CHSH parameters, but thus far have 
found no satisfactory agreement between model calcula- 
tions and measurements. Nevertheless, the density ma- 
trix for the ensemble over the full memory time is poten- 
tially useful for tasks such as entanglement connection, 
as shown by panels C and D where cumulative data are 
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given. Each point at memory time tm gives the violation 
obtained by taking into account all the states generated 
from to Tm- Overall significant violations are obtained, 
namely S+ = 2.21 ± 0.04 > 2 and = 2.24 ± 0.04 > 2 
at Tm ^ 10/is. 

In our experiment, we are able to generate excitation- 
number entangled states between remote locations, which 
are well suited for scaling purposes, and, with real-time 
control, to operate them as if they were effectively polar- 
ization entangled states, which are appropriate for quan- 
tum communication applications, such as quantum cryp- 
tography. Measurements of the suppression ft,*^^^ of two- 
excitation components versus storage time demonstrates 
explicitly the major source that causes the extracted po- 
larization entanglement to decay, emphasizing the critical 
role of multi-excitation events in the experiments aiming 
for a scalable quantum network. The present scheme, 
which constitutes a functional segment of a quantum re- 
peater in terms of quantum state encoding and chan- 
nel control, allows the distribution of entanglement be- 
tween two quantum nodes. But the extension of our work 
to longer chains involving many segments becomes more 
complicated, and out of reach for any current system. 
For long-distance communication, the first quantity to 
improve is the coherence time of the memory. Better 
cancellation of the residual magnetic fields and switching 
to new trap schemes should improve th is p arameter to 
~ 0.1 sec by employing an optical trap [30|, thereby in- 
creasing the rate of preparing the ensembles in the state 
of Eq. Il]) to ~ 100 Hz. The second challenge which 
would immediately appear in an extended chain would be 
the increase of the multi-excitation probability with the 
connection stages. Recently, Ref. [20| has theoretically 
demonstrated the prevention of such growth in a simi- 
lar setup, but its full scalability still requires very high 
retrieval and detection efficiency, and photon-number re- 
solving detectors. These two points clearly show that the 
quest of scalable quantum networks is still a theoretical 
and experimental challenge. The availability of our first 
functional segment opens the way for fruitful investiga- 
tions. 



APPENDIX A: EXPERIMENTAL DETAILS 

Ensembles {LU , LD) are pencil-shaped groups of cold 
Cesium atoms in a magneto-optical trap (MOT) while 
ensembles {RU , RD) are in another MOT, 3 meters 
away. |s), |e)} correspond to the hyperfine levels 

{\6Sy2,F = 4),|65i/2,F = 3),|6P3/2,i^' - 4)}, respec- 
tively. In each MOT, the ensembles U, D are separated 
by 1 mm by way of birefringent beam displacers [27j. 
The MOT is formed at a repetition rate of 40 Hz. In 
each cycle, the MOT is loaded for 18 ms, after which 
the magnetic field is quickly switched off. The trapping 
beams are turned off 3 ms after the magnetic field, while 
the repumping beam stays on for another 100 fis before 



being switched off in order to prepare the atoms in the 
F = 4 ground state \g). 3.4 ms after the magnetic field is 
turned off, trials of the protocol (each consisting of suc- 
cessive write, read, and repumping pulses) are repeated 
with 575 ns period for 3.4 ms. In each trial, the write 
pulse is ~ 30 ns in duration and 10 MHz red-detuned 
from the \g) — > |e) transition. The read pulse and the 
repumping pulse are both derived from the read beam 
(resonant with the |s) — > |e) transition) with 30 ns and 
75 ns duration, respectively. The read pulse is closely fol- 
lowed by the repumping pulse. The read pulse is delayed 
w 400 ns after the write pulse, leaving time for the con- 
trol logic to gate it off, along with the subsequent pulses. 
Independent phase stability measurements show that the 
phase T] drifts in a negligible way, (7r/30) over 500fis cor- 
responding to 870 trials. Some other parameters of the 
experiments are calibrated and listed in table IH 

APPENDIX B: FRINGE VISIBILITY AS A 
FUNCTION OF h^'^^ 

Let us consider that the two pairs of ensembles, U and 
D, have been prepared by heralded detections at Dia, Du 
and Die, Did- Denote by pio, poi, and pu the probabil- 
ity pij to register i photodetection events in field 2m 
and j in field 2^^/ after firing the read pulses. We will 
assume, for simplicity, the various pij are the same for 
both pairs of ensembles. For each of them, the suppres- 
sion of the two-photon events relative to the square of 
the probability for single-photon events is characterized 
by the parameter fe^^^ [211]: 

^ .P1L_ . (Bl) 
PioPoi 

We next relate h'-^^^ to the maximal Cmax and minimal 
Cmin coincidence probabilities between various output 
ports of the detection polarizing beamsplitters {PBS) 
for the left and right nodes at detectors Z?2a , and 
D2c,D2d (see Fig. 1 of the main text). Consider, for 
example, the transmitted ports of the PBS at the L, R 
detectors for the case that the left node has the half-wave 
plate (■|)l set to 0°. In this case, fields 2l(j and 2^^ 
are detected independently, with field 2l]j transmitted 
at the PBS. On one hand, Cmax is obtained for crossed 
polarizers (i.e., {^)r set to 45° at the right node, with 
then field 2im transmitted) and is given to lowest order 
by: 

Cmax^PlOPOl- (B2) 

This term corresponds to the case where only a single 
excitation is distributed in each pair, and each retrieved 
photon is detected from a transmitted port on each side 
L,R. 

On the other hand, the minimum coincidence probabil- 
ity Cmin is obtained for parallel polarizers (i.e., {^)r = 
0° at the right node, with then field 2i^£i transmitted) 
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TABLE I: Noise and Efficiencies 





U 


D 


Field 1 dark count rate 


~ 10 Hz 


~ 10 Hz 


Field 2 dark count rate 


~ 100 Hz 


~ 100 Hz 


Overall retrieval efficiency Pc 


6.4% ± 0.5% 


8.0% ± 0.5% 


Field 2 propagation loss 


68 ± 5% 


68 ± 5% 


Field 2 photon detection efficiency 


50 ± 5% 


50 ± 5% 



and can be written as: 



Pll- 



(B3) 



This term corresponds to coincidences due to photons 
coming from the same pair of ensembles. The smaller is 
the excitation probability, the smaller is this background 
term. 

Taking Eqs. (|B2[) and (jB3p into account, we find that 
the visibility for the number of coincidences as a function 
of the right polarizer angle (i.e., the angle for {^)r) is 
given by: 



V 



c„ 



Cn 



Cn 



1 - 

1 + /l(2) 



(B4) 



Assuming that the visibility is the same in each basis, we 
then find a CHSH parameter S equal to j33j : 



s = 2V2V = 2V2 



I'D 

A minimal value h)^[^^ = 0.17 is thus required to violate 
the CHSH inequality 5 < 2 in the absence of any imper- 
fections except the intrinsic two-photon component. This 
value underlines that this experiment is much more strin- 
gent in terms of minimization of high-order terms than 
previously reported setups. For example, in Ref. [2^ . 
where entanglement between a photon and a stored ex- 
citation is reported, a value of /i*^^-' equal to 0.68 was suf- 
ficient to violate the inequality. The dramatic improve- 
ment reported recently by different groups for the quality 
of the photon pairs emitted by an atomic ensemble was 
thus an enabling step for the practical realization of such 
a more elaborate procedure involving a total of 4 ensem- 
bles reported in the main text. 



1 - fe(2) 

1 + /l(2) 



(B5) 



APPENDIX C: TWO-PHOTON INTERFERENCE 
AND INFERRED OVERLAPS 



two-photon component can lead to coincidences, which 
reduce the visibility. Let us determine the expected visi- 
bility as a function of the two-photon component by way 
of a simple model. 

Consider P„ the probability of finding n photons in 
field 2, and assume the various P„ are the same for both 
ensembles involved. In the ideal case where all ensembles 
have the same properties, the two-photon suppression for 
each field 2 can also be characterized by the same /i'^^) 
parameter used before, which can be written here as: 



2^ 

Pi 



(CI) 



When the half-wave plates {^)l,{^)r are at 0°, the 
fields 2 are detected independently and the probability 
Pmax to register coincidences is given by: 



Pn 



Pi 



(C2) 



When the half-wave plates {^)l,{^)r are rotated to 
22.5°, if the two fields overlap perfectly, the term with 
one photon in each input does not lead to coincidences. 
If we denote by f the overlap, the probability Pmin to 
have one photon in each output is then: 



Pn 



P2 P2 r 
- — + — = 1- 

2 2 ^ 



-/.(2)]^ . (C3) 



From these two probabilities, we find that the visibility 
of the dip in coincidences can be written as: 



Pri 



P-n 



l+^2_/i(2) 



Pr, 



(C4) 



In our case, the measured visibility Vdip is 85 ± 2% for 
the left node and 89±2% for the right one. The measured 
average /i*^^) parameter for this set of data is 0.09 ± 0.01, 
which should lead in the case of perfect overlap to visibil- 
ities Knodoi = 95.5±0.5%. From the measured visibihties 
and this simple model, we can then estimate the overlaps: 
i = 0.89 ± 0.03 for the left node and ^ = 0.93 ± 0.03 for 
the right node. 



For a non-perfect overlap ^ of the field-2 photon 
wavepackets, the visibility of the fringes in the 45° basis 
is decreased by a factor ^"^ . This overlap can be deter- 
mined by two-photon interference, which is implemented 
by mixing the fields 2u and 2d on each side (Right 
and Left) by rotating the half-wave plates (f )l7 (f )i?. by 
22.5°. If the single photon wavepackets are indistinguish- 
able, no coincidences should be observed. However, the 



APPENDIX D: DECOHERENCE TIME OF THE 
STORED EXCITATION 

Residual magnetic fields, which lead to inhomogeneous 
broadening of the ground states levels, is the major lim- 
iting factor of the coherence time Tc of the stored excita- 
tion [l^,!!^!. Consequently, if we neglect dark counts, the 
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conditional retrieval efRciency pc = poi + pio is expected 
to decay exponentially with the storage time tm- 



Pc 



I s 

Pc exp( ) 



(Dl) 



Figure [5] shows an independent measurement of p^ vs. 
tm, with the U and D pairs separated. Fitting the data 
with Eq. (|D1[) gives, for the U and D pairs respectively, 
pI = 7.0% ±0.1% and 8.7% ±0.2%, and Tc = 9.1 ±0.6 /iS 
and 8.5 ± 0.5 /iS. 





FIG. 5: Conditional probability pc of detecting one photon in 
a field 2 for the U (black squares) and D (red circles) pairs, as 
a function of the storage time tm of an excitation. The error 
bars indicate statistical errors. The solid lines are fits using 
Eq. (IdT|| . 



The decay of Pc leads to a similar exponential decay 
of Cij. Cij (i, j = a,b,c,d) are the coincidence count 
rates of two field 2 photons conditioned on two heralding 
field 1 photons defined before. Summing over all Cij used 
in calculating S±, we obtain the total coincidence count 
rates Cs± for the measurement of the Bell parameters S+ 
and S-. Cs±{tm) corresponds to the probability distri- 
bution of the 5'±(tm), and is reflected in the statistical 
error bars /^S± {tm ) ■ The decay of Cs± with tm is shown 
in Fig. [SI Fitting the data with exponential functions: 



s± 



Cs± exp(-TA//T±), 



TM > 0, 



(D2) 



gives T+ — 9.1 ± 0.4 ^s and t_ = 8.1 ± 0.3 ^s, in good 
agreement with Tc. Note that C^^ — 2Cs±{tm = 0), 
since Cs±(r 0) is conditioned on two excitation in 
a same trial, while Cs±(t > 0) is conditioned on two 
excitations created in two different trials: the factor of 
two accounts for the two possible orders of excitations 
('U' then 'D' or 'D' then 'U'). 



APPENDIX E: CONDITIONAL CONTROL AND 
INCREASE IN GENERATION RATE 

As demonstrated in Ref. the conditional control 

of remote memories enables a large enhancement of co- 
incidence rates relative to the case where no logic is im- 
plemented. If the state prepared in one pair of ensembles 



FIG. 6: The tm dependence of the total conditional count 
rates Cs± in the measurement of (a) S+ and (b) S-. The 
horizontal thick lines indicate the size of the memory bin. 
The error bars indicate statistical errors. The solid lines are 
fits using Eq. (|D2)). 



is held up to N trials, the rate for preparing both pairs 
is increased by a factor (2iV + 1) for very low excitation 
probability [25|. Figure [Tl^a) gives the probability pn of 
simultaneously preparing the two pairs. After 17 trials, 
an increase by a factor 34 is obtained experimentally, 
close to the expected value of 35. The gain in the prob- 
ability P1122 of generating the effective entangled state 
is expected to be the same if the coherence time is long 
enough. However, our finite coherence time results in a 
smaller increase of the probability to detect field 2 coin- 
cidences. This increase is given in Fig. [TJ^b), with a com- 
parison to the ideal case of very long coherence time. A 
19-fold enhancement is finally obtained. Let us note that 
the different experimental rates can be obtained from 
these probabilities times the number of trials per second 
(- 2.36 X lOVs). 




FIG. 7: Probabilities of coincidence detection as functions of 
the number of trials A'^ for which the first prepared pair holds 
the state, (a): measured probability pn of preparing the two 
pairs, (b): measured probability P1122 of detecting field 2 co- 
incidences. The green solid line corresponds to the ideal case 
of very long coherence time. Both panels give in addition to 
these probabilities the enhancement factor obtained relative 
to the case without conditional control. 
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APPENDIX F: CORRELATION FUNCTIONS 

^(O",^^), £(45°, ^ii) FOR THE IDEAL EFFECTIVE 
STATE 

In practice, various imperfections lead to deviations 
from the ideal effective state, Eq. (2) in the main text. 
We have developed a detailed model relevant to our ex- 
periment, but consider here only a generic form. Col- 



lective excitations are not shared with equal amplitudes 
between a pair of ensembles because of nonidealities in 
the writing and heralding processes. Likewise, the map- 
ping of atomic states to states of field 2 by the read pulses 
is not ideal. Overall, these various imperfections lead to 
a state \iP2lu,2ru,2ld,2rd) for field 2 given by (neglecting 
multi-photon processes) : 



®(e«D|02,J|l2„^)±e^''°eLD|l2,^)|02„^)) 

= eflc/eflD|02,„>|02,„>|l2H.)|l2«^)±e^''"e"'°ei[;eLD|l2,„)|l2,^)|02Hj|02,,) 

±e*''"eflDeLi/|l2,,)|02,,)|02„,)|l2«^)±e^''"eflC/eLi?|02,,)|l2,^)|l2„„>|02,^) , (Fl) 



where ex is the probability amplitude that a photon is 
created in field 2x ■ The first and second terms in the ex- 
pansion correspond to the cases that the two excitations 
are both retrieved at node "right" and "left", respec- 
tively. Thus the effective state that yields one detection 
event at node "left" and the other at node "right" con- 
sists of the last two terms. After the fields are combined 
by PBSl and PBSi,;, we get the (unnormalized) effective 
state of fields 2l and 2^ 

1^2., 2 Je// = a\H2M±P\V2,H2,) , (F2) 

where a oc e""^ euu^LD and /3 cx e^^" euDiLU- 

From the effective state \'ip2L-2R)eff, we can derive the 
various coincidence probabilities Pij, i,j £ {a, 6, c, d}, 
where {a, b, c, d} refers to the detectors I?2{a,fc.c,ci} for field 
2 in Fig. 1 of the main text. When 9l is fixed at 0°, we 
find (assuming unity detection efficiency) 

Pac = \a\hm^0R 
Pm = |/3|'sin20ji 

Pad = H^COs'Or 
Pba = |/3pCOs20fl 

EiO°,9R) cx Pac + Pbd-Pad-Pbc = -cosi29R) 

(F3) 

irrespective of the ± sign. 

By contrast, when 6*^ is fixed at 45°, we obtain 

Pac = i[l±2|a||/3| cos0cos(9O°-20k) 
+ (|/3|2 - |a|2) sin(90° - 29r)] , 



where cj) — arg{/3)-'arg{a). Let a — cos(p, and /? = simp. 
Denoting S = 45° — 9r, we have 

Pac = ^ |sin2(^| COS0 cos2(5 — cos2(^ sin25] 

Pm = ^^fl ^ |sin2(^| cos0 cos2(5 -I- cos2(/3 sin2(5] 

Pad — ^[-'^ "F |sin2(^| COS0 cos2(5 4- cos2(^ sin25] 

Pbc = ^[1 "F |sin2(^| COS0 cos2(5 — cos2ip sin25] 

E{'15°,9r) cx Pac + Pbd- Pad- Pbc 

= ±\sm2ip\ cose/) cos2(5. (F4) 

From the expression for i?(45°, 9r), we see that the devi- 
ation of \a\ and |/?| from the balanced value, l/\/2, will 
lead to reduction in the visibility of £{45°, 9r) fringes 
and thus the magnitudes of the CHSH parameters S(±). 
We believe that such an imbalance is responsible for the 
results displayed in Fig. 3(b) for E{4:5°,9r) and Fig. 4 
for S'(±) at tm = in the main text, with measurements 
underway to quantify this association. 

Note that another combination of -Py 's given above 
results in 

FiA5°,9R) = -Pac + Pbd+Pad-Pbc 

= cos2cp sin2(5. (F5) 

F{4:5°,9r) allows us to determine cp and thus the magni- 
tude of the coefficients a and /3, independent of(j). Specif- 
ically, the visibility of the i^(45°,6'fl) fringes normalized 
to that of E{0°,9r) fringes yields cos2ip. 
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